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Abstract 

Experimentally, the intermetallic compound Mn4FeGe3 has been recently shown to exhibit en- 
hanced magnetic properties and spin polarization compared to the MnsGea parent compound. 
The present ab-initio study focusses on the effect of Fe substitution on the electronic and magnetic 
properties of the compound. Our calculations reveal that the changes on the Fermi surface of 
the doped compound are remarkable and provide explanations for the enhanced spin-polarization 
observed. Finally, we show that it is indeed possible to tune the degree of spin-polarization upon 
Fe doping, thus making the Mni_a;Fej;Ge3 intermetallic alloy very promising for future spintronic 
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PACS numbers: 71.15.Mb;71.20.Lp; 72.25.Ba; 72.25.Hg 



The ferromagnetic intermetallic MnsGes compound is a very promising spin injector: 
it can be easily grown epitaxially on the Ge(lll) substrate Q, Q, and can provide spin- 
polarization of about 42 % while preserving ferromagnetic ordering up to a Curie tempera- 
ture (Tc) as high as 296 K. Very recently, T. Y. Chen et al. have shown that replacing one 
Mn atom/f.u. with Fe may lead to a remarkable effect on the spin-injection properties. [4] 
This immediately opened the possibility for Spin Polarization Engineering (SPE) of MusGes 
by means of Fe doping. The reasons for the remarkable spin polarization enhancement in 
Mn4FeGe3 are not clear yet and await theoretical explanations. 

Density Functional calculations were performed using the VASP package within the Gen- 
eralized Gradient Approximation (PBE-GGA). 8] PAW pseudopotentialsj^, Q were used 
for both Ge and TM-atoms: the semicore 3p states are considered as valence (core) for Mn 
(Fe); the Sci states are frozen in the core for Ge. The kinetic energy cutoff used for the wave 
functions was fixed to 350 eV. (4,4,6) P-centered fc-points were used for the self-consistent 
cycle, while a (12,12,14) k-point grid was used for the calculation of the Fermi velocity. All 
the atomic internal positions as well as the volume and shape of the unit cell were relaxed 
minimizing the ab-initio stress and forces. 

The Mn4FeGe3 compound can be represented as a solid solution of Fe in MnsGes, i.e. 
Mn5_j.Fe^Ge3 with x = 1. For x = 0, at ambient conditions, MusGcs crystallizes in the 
hexagonal DSs type (space group mem) with a unit cell containing 16 atoms: 10 Mn 
atoms in two inequivalent sites identified using the Wyckoff notation [4M/ in 4((i); 6M// 
in Q{g)] and 6 Ge atoms in 6{g) sites. The description of the crystal structure has been 
given elsewhere, jul For x = 5, FesGes is also hexagonal but belongs to the D82 symmetry 
type (space group P63 /mme), so that the mutual solubility of the end members could be 
limited by their different crystal structure. The limit of solubility appears to occur at 
x=l, or slightly beyond, [sl resulting in the Mn4FeGe3 compound that still preserves the 
same crystal structure as Mn^Ge^. Therefore, Mn4FeGe3 can be described using an ordered 
supercell containing 16 atoms, i.e. 8 Mn, 2 Fe and 6 Ge atoms. 

We start our study considering one Fe atom at both the M/ or M// site, and we found 
that Fe prefers to occupy the M/ site with an energy gain of 280 meV/cell with respect to 
the Mil site: thus, in agrement with experiments, si, m, B] the smaller and less electroposi- 
tive iron atom substitutes preferentially Mn on the 4{d) sites. |l2|| This not unexpected since 



occupation of the Mj site allows the metallic atoms to be closer 



consistently with the Fe 



TABLE I: Structural data (lattice constant a and c/a ratio, magnetic moments and heat of forma- 
tion) for the structures considered compared with experiments where available. 



System a(A) c/a (fJ-B) HMu (f^B) ^J'Ge (^J-b) fJ^T ()UB/f.u.) AHf{eV/atom) 

MnsGes 7.142 0.697 2.22 3.11 -0.16 2.70 -0.140 

Exp[13] 7.184 0.703 1.96 3.23 2.60 

Mn4FeGe3 7.134 0.690 1.75 (Fe), 2.35 (Mn) 3.05 -0.15 2.59 -0.145 

Exp[4] 7.184 0.696 2.35 

Exp[7] 7.138 0.702 1.55 2.45 2.10 

FesGes 6.967 0.692 1.58 2.17 -0.11 1.85 -0.061 



atoms having smaller atomic radius than Mn. The 4{d) atomic positions (in internal coor- 
dinates) lie in two different planes along the hexagonal c-axis: (1/3,2/3,1/2), (2/3,1/3,1/2), 
(1/3,2/3,0), (2/3,1/3,0) and they can be occupied by two iron atoms in three different and 
non-equivalent configurations. We find that, in the lowest-energy state, the Fe atoms are 
located on the sites belonging to the z = or z = c/2 symmetry equivalent planes. 

In table [H we show the calculated lattice constants and magnetic moments of MnsGea, 
Mn4FeGe3 and the hypothetical FesGes compound in the DSs crystal structure. The for- 
mation energy {AHj) is defined with respect to reservoirs of atoms in the pure bulk phases 
7-Mn AFMl, FM bcc-Fe and diamond Ge; we have chosen the sign in such a way that stable 
compounds have negative formation energies. For Mn^Ge^, the theoretical lattice constants 
a and c are slightly underestimated with respect to the experimental values (0.5 % and 1.4 
%, respectively). The calculated magnetic moments correctly reproduce the relative magni- 
tude of the M/ and Mjj moments, with the former carrying smaller moments, in agreement 
with experiments and previous theoretical calculations. U, Q, 1^ A small induced anti- 
ferromagnetic polarization is present in the cell, mostly localized on Ge sites. 11, Q] 
This is true for all three cases considered. Upon doping, the volume of the unit cell de- 
creases by 1.4 %; the c-axis shrinks faster than the a-axis, thus decreasing the c/a ratio. 
The M/-M7 separation along the 2;-axis {i.e. c/2) affects the length of the c-axis most: the 
smaller atomic size of Fe compared to Mn (the Fe and Mn covalent radii are 1.25 and 1.39 A, 
respectively [12:]) is clearly responsible for the contraction of the c-axis upon Fe substitution. 
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Our calculated lattice constants agree well with experimental values. [4, tI] The magnetic 



moments at the M/ sites are distributed as follows: Fe carries the lowest magnetic moment 
(1.75 fis), while Mn shows a slightly increased moment compared to pure MnsGes (2.35 
fis)- On the other hand, the M// magnetic moments are only marginally affected by Fe 
substitution. The experimental values shown in Tab. [T] are obviously average values for the 
two independent sublattices: the slight discrepancy between experimental and calculated 
values could be ascribed to disorder effects which are not taken into account in the present 
calculation. 

It is interesting to explore the x = 5 doping limit in the D8s phase. In the FesGca phase, 
there is a significant decrease of the length of both the a— and c— axis by 2.4 and 3.1 % 
respectively, while the c/a ratio is almost unchanged compared to the x = case. The 
iron magnetic moments decrease to 1.58 and 2.17 fis at M/ and Mjj sublattice, leading 
to a sizable reduction of the total moment in the cell. Finally, concerning the stability of 
the compounds, we note that MnsGes and Mn4FeGe3 are quite stable (Aif/ are -0.140 and 
-0.145 eV/atom respectively), whereas FesGcs has only a slightly negative formation energy 
(-0.061 eV/atom). 

We now focus on the electronic structure of the intermetallic compounds considered. 

n 

According to Ref. |6|, strong metal-germanium and metal-metal covalent bonding is present 
in this phase. Broad, primarily sp bonding and antibonding bands form due to covalent 
interactions between Mn and Ge. Since there are three germanium atoms per formula unit, 
each contributing four orbitals with a spin degeneracy of two per orbital, these bands can 
accomodate 24 electrons. Assuming formal valences (Mn,Fe)"''^ at 4((i), (Mn,Fe)"^^ at 6(5'), 
Ge^~ at 6(5^),^] the bonding band is filled with six electrons from the two 4{d) metal atoms, 
six from the three 6{g) metal atoms, and twelve from the three metalloid atoms (in total 
24). At the Fermi level, a large density of d-like states is present. Therefore, we expect 
that there will be sheets at the Fermi surface with large and "heavy" -like areas, made 
up predominantly of electrons relatively localized on partially occupied d-\ike states and, 
thus, with large effective masses and low Fermi velocity. On the other hand, there will be 
also small and "light" -like areas of the Fermi-surface sheets with sp character, made up by 
highly mobile electrons with delocalized wave functions, small effective masses, and high 
Fermi velocity. Thus, it is expected that the d-likes state dominate the density of states at 
the Fermi level, N(Ej), while sp electrons do contribute substantially to the average value 
of the Fermi velocity. 
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FIG. 1: (color online) Electronic density of states (DOS) of MnsGea, Mn4FeGe3 and FesGes. Total 
and projected DOS onto metal 4((i) and 6{g) sites are shown in top and bottom panel. The sp-DOS 
has been multiplied by a factor of 4 for clarity. 

Fig. [1] shows the spin resolved total DOS (upper panels) and the DOS projected onto 
(i-states of 4(d) and 6(g) metal atoms as well as onto sp states (lower panels). From left to 
right, the iron concentration varies from x=0,l to x=5. Let us focus on the region around the 
Fermi level for MusGca (upper panel): while the spin up DOS has an almost structureless 
plateau (there is a small valley in the DOS curve just above the Fermi level), the spin down 
DOS has much more fine structure around and immediately above the Fermi energy. From 
the lower panel, we see: i) the sp electrons are mainly located in the region below -2 eV 
(bonding states) and above 1.5 eV (antibonding states), supporting the picture given above; 
ii) N(Ej) is dominated by d states rather than by sp states, in both spin channels, suggesting 
that the d electrons are mainly responsible for the spin polarization; iii) the number of sp- 
electrons with spin-up, Njp, is significantly larger than those with spin-down, N^^, suggesting 
that the average value of the spin-up Fermi velocity is larger than the corresponding spin- 
down value. (iTj] In the majority component, the d-A(d) DOS is larger than the d-6{g) one, 
while they are practically equal in the minority component. Upon Fe doping (electron-like 
carrier doping), the Fermi level shifts upwards to accomodate the extra electron per formula 
unit and now lies at the bottom of the small valley of the majority DOS and on top of 
a small peak in the minority DOS: as a result, N^^ decreases, and N^- increases, leading to 
a more negative spin-polarization Pq, which is proportional to N^^-N-^ (see below). These 



TABLE II: Spin polarizations for different x = 0, 1, 5. See text for further details. 
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Mn4FeGe3(x = 1) 
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modifications are mainly driven by the d-4{d) states and are consistent with a simple rigid 
band picture. On the other hand, there are negligible changes in Nj^ and N^^ at Ej. 

In the high doping limit, we observe that the majority DOS is hardly modified compared 
to the lower Fe content case, while the minority component becomes more structured; in 
fact, the Fermi level is further shifted upwards: it still lies at the bottom of a valley in the 
majority DOS, while it is pinned at a high peak in the minority component, leading to a 
very large value of the total DOS at Ej. The large Fe content mainly affects the minority 
spin DOS and strongly modifies the overall transport properties through large modifications 
of the minority spin component. As before, the sp-DOS at Ep is marginally affected. 

Let us now consider the spin-polarization (SP) of the different compounds. A common 
definition for SP is: jisl 

where N and vp are the density of states (DOS) and average Fermi velocity of electrons 
with spin-up, spin-down at Ef, respectively. Note that Eq. [T] can be written as follows: 
Pn = {xoVq - l)/(a;ol/o + 1) with xo = N\Ef)/N^{Ef) and = 4/^. Due to the 
hexagonal symmetry, the Fermi velocity can be decomposed into components which are 

T I T I 

parallel and perpendicular to the basal plane, Vp]^^ and Vp^ respectively. The same is true 
for yl. 

For n = 0, the SP calculated using Eq. [T] corresponds to spin-resolved photoemission 
measurements, while higher orders of P„ correspond to SP as measured in transport exper- 
iments, such as PCAR (point contact Andreev reflection ) and TJ (tunnel junction), in the 
ballistic (Pi) or diffusive regime (^2)- 161] In Tab. [Tll we show the calculated xq and yo at 
Ef for the three compounds, and the corresponding P„ for ?t,=0,1,2. Upon iron doping, we 
see that: i) <N^ {i.e. xq < 1) in all compounds considered and, in particular, <CN-'- 
for X = 5. ii) For both the parallel and perpendicular component, v'p > Vp {i.e. yo > 1), 



except for FesGea where the opposite is true. The average velocity for spin up Vp is larger 
than Vp, because the conductivity is dominated by the sp electrons, and in all Mn5_a.Fea;Ge 
compounds N|p >N^p. FesGes is special, since N;^ is exceptionally large (see Fig. 1). 

Clearly, point i) leads to negative Pq, with the absolute values increasing with the Fe 
content; i.e. at the Fermi-level, the spin up density is smaller than the spin down density, 
and this behaviour is enhanced by Fe doping. However Pi is greater than Pq, and P2 is 
greater than Pi, except for FesGes. This is related to ii), i.e. the spin up velocity is 
exceeds the spin down velocity {i/q > 1) over-compensating the lower spin-up versus spin- 
down density at Fermi- level. The exception is FesGea where i/q < 1, and hence P slightly 
decreases with n { Pq > Pi > P2). 

In summary, we have studied the hexagonal phase (DSg) of the Mus.^Fe^^Ge intermetallic 
alloy, focusing on the role of the Fe substitution on the structural, electronic and magnetic 
properties of the compound. We found that: i) Fe substitution preferentially occurs at the 
4((i) site of the DSg structure; ii) substitution of one Mn atom/f.u. in MusGes enhances the 
spin polarization Pq in agreement with experiment, but iii) substitution also decreases Pi 
and P2. The last observation is at variance with the experimental results. j4] The reason for 
this disagreement deserves further studies that should include disorder as well as interface 
effects, which are beyond the purpose of the present study. On the other hand, our ab-initio 
calculations indeed confirm that is possible to tune the degree of spin-polarization in the 
different transport regimes upon Fe doping. Remarkably, an even larger spin polarization 
than that of Mn4FeGe3 can be achieved by further increasing the Fe content, thus making 
the Mni_a;Fea;Ge3 intermetallic alloy very promising for future spintronic applications. 
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